function [max_ecg min_ecg sd_ecg] = heartrate_hist(hrc_data)

    HR_modified = hrc_data(3:length(HR_real)-3);
    max_ecg = max(HR_modified);
    min_ecg = min(HR_modified);
    sd_ecg = std(HR_modified);
    plot(HR_real);
    axis 'tight'
